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ABSTRACT 

We present the results of our follow-up observation program of 7-ray sources detected by the 
Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. 26 blazars and 39 
sources unidentified at other wavelengths were targeted at IRSF 1.4 m telescope equipped with 
the SIRIUS/SIRPOL imager and polarimeter. i7-band magnitudes of the blazars at the epoch 
of 2010 Dec - 2011 Feb are presented, which reveal clear flux variation since the Two Micron All 
Sky Survey observations and can be useful data for variation analyses of these objects in longer 
periods. We also find that nearly half of the 7-ray blazars are highly (>10 %) polarized in near- 
infrared wavelengths. Combining the polarization and variation properties, most (~90 %) of the 
blazars are clearly distinguished from all other types of objects at high Galactic latitudes. On the 
other hand, we find only one highly polarized and/or variable object in the fields of unidentified 
sources. This object is a counterpart of the optical variable source PQV1 J131553. 00— 073302.0 
and the radio source NVSS J131552— 073301, and is a promising candidate of new 7-ray blazars. 
From the measured polarization and variation statistics, we conclude that most of the Fermi/hAT 
unidentified sources are not likely similar types of objects to the known 7-ray blazars. 



Subject headings: BL Lacertae objects: general - 
jets — polarization — quasars: general 

1. Introduction 

A new era of blazar studies has arrived with 
the advent of the Fermi Gamma-ray Space Tele- 
scope (Fermi). Fermi has been carrying out an all- 
sky survey with its main instrument Large Area 
Telescope (LAT) since the science mission phase 
started in 2008. The fir st Fermi-LAT catalog 
(1FGL; [Xbdo et alJlioiOah lists 1451 7-ray detec- 
tions, in which 821 sources are associated (or iden- 
tified) with objects found at other wavelengths. 
The majority of the associated objects are ac- 
tive galactic nuclei (AGNs) dominated by blazars, 
while other extragalactic sources as well as Galac- 
tic objects such as pulsars and supernova remnants 
make smaller contributions. The remaining 630 
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gamma rays: galaxies — galaxies: active — galaxies: 

sources are left unidentified in 1FGL. 

Blazars are considered to be AGNs whose jets 
are aligned with the observer's line-of-sight. Emis- 
sion from the relativistically boosted jets domi- 
nates observed flux, resulting in two broad peaks 
in the spectral energy distribution (SED); one at 
radio to X-rays arising from synchrotron emission 
of accelerated, high-energy particles, and another 
at X-ray to 7-ray arising from the inverse Comp- 
ton scattering of the lower energy photons. Be- 
cause of these emission mechanisms, blazars are 
characterized by strong radio, X-ray, and 7-ray 
radiations as well as high polarization and varia- 
tion across the entire SED. Since blazars dominate 
the 7-ray sky at high Galactic latitudes, the Fermi 
all-sky survey is expected to shed new light on this 
relatively rare and poorly understood population. 

An investigation of blazar radiation mecha- 
nisms (hence the intrinsic SED) is important not 
only for revealing the nature of blazars them- 
selves, but also for measuring the extragalactic 
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background light (EBL). Since high-energy pho- 
tons from blazars interact with optical to near- 
infrared (IR) EBL in the intergalactic space, 
observations of distant blazars can be used to 
infer the EBL spectrum if the intrinsic blazar 
SED is precisely known. With this indirect 



method, lAharonian et al.l (|2006j ) obtained the 



significantly-lower upper limits of near-IR EBL 
than those d erived from the direct measure- 
ments by, e.g., Matsumoto et al. ( 20051 ). Recently 
Matsuoka et al.l (|201ll ) succeeded in a direct mea- 
surement of optical EBL by re-analyzing the Pi- 
oneer 10/11 data, which is on the smooth ex- 
tension of the n e ar-IR upper limits obtained by 
Aharonian et al.l (|2006h . 

In this paper, we present the results of our near- 
IR follow-up observation program of 7-ray blazars 
and unidentified sources in 1FGL. Despite of the 
wealth of information there, near-IR wavelength 
of the whole AGN population is sti l l poorly un - 
derstood (e.g.. iMatsuoka et al.ll2007l I2008ll2012t) . 
We aim to quantify the most distinct features of 
blazars in near-IR wavelengths, i.e., polarization 
and variation, in part as a benchmark for future 
observations. At the same time, we explore the 
nature of 1FGL unidentified sources by compar- 
ing their near-IR properties to the known 7-ray 
blazars. Dominance of blazars in the 7-ray sky 
implies that some of the unidentified sources at 
high Galactic latitudes are similar objects missed 
in the past surveys at other wavelengths. Reveal- 
ing the origin of these unidentified 7-ray emissions 
is of greatest importance, therefore many follow- 
up studies are being dedi cated to this subject (e.g., 
Ackermann et al]l2011bl ) . We aim to provide com- 
plementary information to these follow-up pro- 
grams from a near-IR view point. 

While the seco nd Fermi-LAT catalog (2FGL; 
Nolan et al.ll201 21 ) has already been released, this 
paper is based on IFGL for consistency with the 
sample selection of presented observations. We 
discuss the new 2FGL identifications of the sam- 
ple later. Magnitudes are presented on the Vega- 
based system throughout this paper. 

2. Observations and Reduction 

The sample consists of blazars and unidentified 
sources extracted from IFGL. We selected targets 
at high Galactic latitudes (|6| > 20°), with small 



position errors (< 5 arcmin with 95-% confidence 
in most cases), and with bright Two Micron All 
Sky Survey (2MASS3) magnitudes (# 2 mass < 15 
mag) for the blazars whose accurate coordinates 
are known. The objects listed in the Fermi -LAT 
first AGN catalog QLAC; lAbdo et al.ll2010bl) were 
removed from the sample of unidentified sources. 
The polarimctric observations were performed 
with the Infrared Survey Facility (IRSF) 1.4-m 
telescope at the South Africa Astronomical Ob- 
servatory, Sutherland. We used the near-IR imag- 
ing camera SIRIUS (Sim ultaneous Infrared Im- 
ager for Unbiase d Surv ey; iNagashima et al. 19991 
Nagavama et al.l 12003} ) which is equipped with 



three 1024 x 1024 HgCdTe arrays with the pixel 
scale of 0".453. It was designed to obtain 7'. 7 x 
7'. 7 images in the three bands, J (1.25 ^m), H 
(1.63 /jm) and K s (2.14 /im), simultaneously. The 
instrument provides polarimctric capability called 
SIRPOL, with an achromatic (1-2.5 /im) wave- 
plate rotator unit and a polarizer. Source fluxes 
(-Fooo, -F225, -F450, ^675) in the four wave-plate 
angles, 0°.0, 22°. 5, 45°.0, 67°. 5, are measured in 
consecutive exposures. The uncertainty of a mea- 
sured polar ization degree is est imated to be less 
than 0.3% (|Kandori et al.ll2006l ). 

We summarize the observed targets in Table 
Q] (for blazars) and [2] (for unidentified sources). 
The observations were carried out in 2010 De- 
cember and 2011 February. Mean total exposure 
times were 30 min for a blazar and 60 min for an 
unidentified source, divided into single exposures 
of 20 sec between which the telescope pointing was 
dithered by 20 arcsec. Typical seeing during the 
observations was ~1" .4. Dark and twilight-flat im- 
ages were obtained before and/or after each night 
of observations. 

Data reduction was performed in a standard 
manner with a dedicated package SIRPOL of the 
Image Reduction and Analysis Facility (IRAH3), 
including dark subtraction, flat fielding, and sub- 



1 This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the Uni- 
versity of Massachusetts and the Infrared Processing and 
Analysis Center / California Institute of Technology, funded 
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the National Science Foundation. 
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versities for Research in Astronomy (AURA) under coop- 
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Table 1: Observation journal of blazars 
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"The observations were carried out from 2010 December to 2011 February. 

limiting magnitudes below which photometry error in a single wave-plate image is less than 0.05 mag. 
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Table 2: Observation journal of unidentified 
sources 
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"The observations were carried out from 2010 December to 

2011 February. 
b Limiting magnitudes below which photometry error 

in a single wave-plate image is less than 0.05 mag. 



stituting bad pixels. Photometric calibration was 
achieved by referring to several 2MASS sources 
within each observed field. While simultaneous J, 
H, and K s band images were obtained, we use H- 
band images in this study because of their highest 
quality (signal-to-noise ratios). We show an exam- 
ple of a reduced P-band image in the left panel of 
Figure Q] 



We used the Source Extractor (Bertin fe Arnouts 



1996), version 2.5, for aperture photometry of de- 
tected sources. Aperture sizes were determined as 
twice the mean full widths at half maximum of 
stellar profiles in each image. We required pho- 
tometry errors to be less than 0.05 mag for polar- 
ization measurements, which defines our limiting 
magnitudes P]gg P listed in Table [T] and H The 
average values are PiRgp = 15.4 and 16.0 mag 
for the fields of blazars and unidentified sources, 
respectively. Polarization degree P ra w was then 
calculated from measured fluxes Fooo, P225, P450, 
F 6 75 as follows: 

Q = -Fboo — -F450 ) 

U = F225 — -^675 > 
I = (-PoOO + -F450 + F225 + ^675)/2, 



Pr: 



100 x 



VQ 2 + u 2 



Here Q and U are the Stokes parameters and I is 
total intensity. A source with strong polarization 
stands out in Q/I or U/I images as demonstrated 
in the right panel of Figure [TJ 

Since Fooo, F225, F450, ^675 are not observed 
simultaneously, temporal variation of the atmo- 
spheric condition causes additional polarimetry er- 
rors to those from the standard sky background 
noise. Hence we estimated polarization error AP 
in a given field from the median P raw value of 
all the detected sources, assuming that most of 
those sources at high Galactic latitudes have no 
intrinsic polarization. Then true polarization de- 
grees were derived by de- biasing measured P raw 
(jWardle fc Kronberdll97i : bllowing 



^W^aw-AP 2 . 

Thus the presented values of polarization and its 
error should be regarded as the lower and upper 
limits, respectively, considering the above assump- 
tion on no intrinsic polarization for most of the 
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.(a) Intensity 



(b) Stokes U/l 



* 



Fig. 1. — IRSF/SIRPOL .ff-band image of intensity (panel a) and Stokes U/I (b) of a field around the blazar 
1FGL J0538.8-4404. The blazar is marked with the boxes. 



detected sources. We provide upper limits of po- 
larization for objects with P law < AP or P < AP. 

3. Results and Discussion 

We present the results of photometry and po- 
larimetry measurements for the blazars in Table [T] 
Their polarization and variation since 2MASS ob- 
servations, as well as those for all other sources de- 
tected in the blazar fields, are plotted in Figure [2] 
Strong variability of the blazars is evident; 22 out 
of 26 blazars (85 %) have |#irsf-#2MASs| > 0.25 
mag while only one of the other detected sources 
shows such variation (~0.3 mag; it is likely a con- 
taminating normal star). The blazars are also 
characterized by pronounced polarization; 10 out 
of 25 blazars (40 %; polarization was not mea- 
sured for a blazar associated to J0038.4— 2504 due 
to a large (> 0.05 mag) photometry error) have 
P > 10 %, while all other types of objects are 
much less polarized. The presented fraction, 40 
%, represents a useful benchmark for planning 
future optical/near-IR polarization follow-up pro- 
grams of 7-ray selected blazars. 

From the above results, we can derive an ex- 
pected number of blazars which should be discov- 
ered by their high polarization and/or variation in 
the fields of unidentified sources if they are similar 
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Fig. 2. — Polarization and variation of 326 sources 
detected in the blazar fields. The diamonds repre- 
sent the blazars associated to 1FGL sources while 
the dots represent all other objects. The arrows 
denote upper limits of polarization. 
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but unknown blazars. In total, 23 out of 26 blazars 
(88 %) have pronounced polarization (P > 10 %) 
or variation (|i/iRSF—-ff2MASs| > 0.25 mag). Since 
39 unidentified sources were observed, we would 
have ~35 objects with such distinguishable prop- 
erties under the above assumption. However, ^35 
% of them would be fainter than our limiting mag- 
nitude H^ SF ~ 16.0 mag based on the 2MASS 
magnitude distribution of 1FGL blazars. Further- 
more, ~32 % of them would be outside SIRPOL 
field-of-view (7 '.7 x 7 '.7 arcmin 2 ) when 7-ray po- 
sitions are used as the telescope pointing center, 
based on the distribution of distance between 7- 
ray positions and associated blazars. Considering 
these restrictions, we expect to find 15 blazars out 
of 39 fields of unidentified sources if their 7-ray is 
indeed emitted by a similar population to 1FGL 
blazars. 

Figure [3] shows polarization and variation of 
608 objects detected in all the fields of uniden- 
tified sources. While three objects are found to 
have high variation (|i?iRSF ~ #2MASs| > 0.25 
mag), two of them with P < 7 % and |i?rRSF — 
#2MASs| ~ 0.3 mag are likely contaminations; it 
is consistent with one contamination out of 326 
objects in the blazar fields (Figure [5]). On the 
other hand, the object at P ~ 9 % and |i?iRSF — 
#2MASs| ~ 1-2 mag is a promising candidate of 
new 7-ray blazars. It is found in J1315.6— 0729 
field, at R.A. 13 h 15 m 52 s .98, Dec. -07 d 33 m 01 s .99 
(J2000.0) with H-band brightness H mSF = 14.1 
mag and H2MASS = 15-3 mag. Its counterparts 
are found in the NEeH; the opt ical variable source 
PQV1 J131553.00-073302.0 (|Bauer et all l200flh 



and the radio source NVSS J131552-073301. We 
plan to carry out a spectroscopic follow-up obser- 
vation of this object in the near future. 

Except for the above blazar candidate, we 
found no clear sign of variable or polarized ob- 
jects in the fields of unidentified sources. The ap- 
parent inconsistency between the expected num- 
ber of blazar candidates as estimated above (15) 
and the actual number (1) indicates that 7-ray 
emission of the unidentified sources arises from 
other types of objects than known 1FGL blazars. 
They could be non-blazar active galaxies, star- 



3 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronau- 
tics and Space Administration. 



burst galaxies, or Galactic objects such as pul- 
sars and supernova remnants at relatively high 
Galactic latitudes, as well as active galaxies with 
jets but without strong polarization and variabil- 
ity at near-IR wavelengths. In this regard, it is 
noteworthy that some of them are already iden- 
tified (or associated) in the latest 2FGL catalog 
and the Fermi-LAT seco nd AGN catalog (2LAC; 
lAckermann et al ] l2011ah as summarized in Table 
[3J Many of them are pulsars, which is consistent 
with our observation results. While J2330.3— 4745 
is associated to the blazar PKS 2326—477, their 
separation is relatively large and the blazar is out- 
side the field-of-view of our SIRPOL observation. 

We are grateful to the IRSF team at Nagoya 
University, Kyoto University, and National As- 
tronomical Observatory of Japan for great assis- 
tance provided during the observations. This work 
was supported by Grant-in-Aid for Young Scien- 
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